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Abstract 



Solutions to gravity with quadratic Lagrangians are found for the 
simple case where the only nonconstant metric component is the lapse 
N and the Riemann tensor takes the form R^j = —kikj, i,j = 1, 2, 3; 
thus these solutions depend on cross terms in the Riemann tensor and 
therefore complement the linearized theory where it is the derivatives 
of the Riemann tensor that matter. The relationship of this metric to 
the null gravitational radiation metric of Peres is given. Gravitaional 
energy Poynting vectors are construcetd for the solutions and one of 
these, based on the Lanczos tensor, supports the indication in the 
linearized theory that nonnull gravitational radiation can occur. 



1 Introduction 

In general relativity it is thought that the assumption of the field equations 
is sufficient to ensure that gravitational radiation is null (i.e. travels at the 
speed of light). It is not necessarily the case that gravitational radiation is 
null but in general relativity four facts support this view. The first is 
that in the linear approximation, with the harmonic gauge, the field equa- 
tions reduce to Oh a b = 0, which has only null wave solutions. The weak 
field linear approximation may no do justice to the non;inearity of the fielf 
equations; for example, on a de Sitter background, with cosmological con- 
stant A, perturbations obey the Fierz-Pauli equation (□ — m 2 )h a i, = 0, with 
mass m 2 = 2A/3 []. The second fact is that exact solutions of Einstein's 
field equations which are good models of gravitational radiation represent 
null radiation. The third is that the charateristics (or shock wave soltuions) 
of the field equations are null; however, this is not a particularly compelling 
reason - for example, the characteristics of the massive Klein-gordon equa- 
tion are null Q. The fourth is that as a consequence of choosing a metric 
geometry (where V c G a b = 0), there is the equation D^ab = 0, which appears 
to be massless; this has led to the construction of nonmetric massive theories 
of gravity obeying the equation (□ + m)g a b = 0. 

One of the most important variants of general relativity is the theory 
with quadratic terms added to the Lagrangian. Einstein's field equations 
of general relativity can be derived from Hilbert's Lagrangian 0]. This La- 
grangian was generalized to include terms of higher order in the early days 
of relativity, by Pauli || , Bach || and others. Recently @ H , the relevance 
of quadratic Lagrangian theories to the shape of galaxies has been discused. 



1 footnote added 1999, this calculation is done section 6 gr-qc/9812091 
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When quadratic Lagrangian theories are linearized 0] [ 10 1, it is found that 
the field equations reduce to two wave equations with masses fixed by the 
values of the coupling constants. In one case the massis positive, leading 
to slower-than-light waves, and in the other it is negative, leading ot faster- 
than-light waves. The initial value formulation of quadratic Lagrangian 



theory [11], suggests that these masses are present in the full nonlinear the- 
ory. Here the existence of exact nonnull wavelike solutions to the quadratic 
Lagrangian theory is investigated. The method used is to start with a mod- 



ification of the Peres' wave [12], where the modification entails changing 
from Peres' dependence on a null coordinate to a nonnull coordinate. Using 
of coordinate transformations, when the coordinate is timelike the metric 
is shown to be equivalent to a metric with the lapse N as the only varying 
component of the metric, and when the metric is spacelike the metric is 
shown to be static. The Peres metric itself has been discussed in the con- 



text of quadratic Lagrangian theories by Buchdahl |T3|] and Madsen [ 14| ; the 
metric they use has cylindrical symmetry and depends on a null coordinat, 
whereas the solutions here do not have these properties. The conventions 



used are those of Hawking and Ellis [ 15 1 . 



2 Quadratic Lagrangians 

The action is taken to be of the form 

S = J a \ [£ m + k~ 2 {R - 2 A) + bC 2 + P R 2 ]dx 4 



(2.1) 



where C m is the matter Lagrangian, the R term is the Hilbert Q action, 
A is the cosmological constant, the C 2 term is the Bach || action with 
C 2 = C a bcdC abcd being the square of the Weyl tensor, and the R 2 term is the 
Pauli H action. Such nomenclature is historically simplistic; for example, 
Pauli was considering theories with nonmetric connection. Linearization j|] 
L0| shows that the coupling constants can be identified with the masses 
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(2.2) 



For A = 0, there are eight physical degrees of freedom, two for the massless 
graviton corresponding to the Hilber term, one massive scalar corresponding 
to the Pauli term, and five massive spin 2 poltergeists corresponding to the 
Bach term. The sign of these terms can also be considered from the point 



of view of cosmological stability 16] 
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Varying 2.1 gives the field equations 

K-\G ab + Ag ab ) + M ab + X ab = T c 



ab i 



where 



Cafe — Rab — -gabR, 



(2.3) 



(2.4) 



M ab = -2baR ab + 2 p + - « ;a6 - 2p - - ffaftD-R, (2.5) 



X ab = 4W? ocrf6 i? cd - 2 I p 



26 



RR ab 
2b" 



R< 



The trace of 2.3 



is 



k~ 2 (4A- R) -6pnR = T, 



(2.6) 



(2.7) 



with T = T c c , which leads to the alternative form of the field equations 

K 2 {Rab — 9abA) + X ab + Z ab = S ab , (2.8) 

where 

Z ab = -2bOR ab + (p + t\ [ 2 R. ab + ^Di?], (2.9) 



(2.10) 



and it is this form of the field equations that is used here. The stress tensor 
is taken to vanish or to be that of a perfect fluid, 



Sab = (jj>+ v)U a U b + -(a* - v)g a b, 



(2.11) 



where p and v are respectively the density and the pressure of the fluid. U a 
is a unit timelike vector field, 



U a U a = -1, 

with acceleration, expansion and projection, 

A a = U a b u b = U a , 



(2.12) 





■ab 



(2.13) 



g ab + Jjajjb^ 
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respectively. The first and second conservation equations are 

fi a U a + fa + v)Q = 0, 

(ji + u)A a + i ab u b = 0, (2.14) 

respectively. 

3 Poynting Vectors 

In a nonvacuum space-time the speed of a gravitational wave can be mea- 
sured relative to the matter present. For example, in a cosmological model 
it can be compared to the speed of the comoving fluid velocity vector; on a 
galactic scale the speed of a gravitational wave emitted from the center of 
the galaxy could be estimated be comparing it to the speed of co-occurring 
electromagnetic radiation, or possibly by direct methods. The metrics that 



will be here have the property O, where the first conservation implies that 
the density of a fluid is at rest in this coordinate system, so that there is 
a possiblity of measuring the speed of a gravitational wave with respect to 
this fluid. So far the field equations coupled to a perfect fluid have proved 
intractable, and only solutions with zero stress have been found; thus for a 
timelike wave it is possible to move to a coordinate system where the wave 
appears to be at rest, and as there is no matter present this criterion cannot 
be used to measure wave speed. Another criterion is needed and the one 
suggested here is to construct Poynting vectors representing the speed of 
energy transfer: whether the Poynting vector is timelike, null or spacelike 
will be taken to support the belief that the gravitational wave is timelike, 
null or spacelike. If the principle of equivalence holds, then the energy of the 
gravitaional field appears to be nonexistent to an observer in free fall, and 
consequently a Poynting vector for gravitaional energy gives at best only an 
indication of the real energetics involved. To sum up: whether a gravita- 
tional enery Poynting vector is timelike, null or spacelike can be regarded 
as measuring whether the gravitational wave is timelike, null or spacelike; 
alternatively it can be considered as a measure of the appropriatenees, or 
otherwise, of the Poynting vector. 

There are tensors constructed out of products and derivatives of the 
Riemann tensor, which can represent the square of gravitational energy |L7]] 
There are a large number of possible combinations pll. An often-used com- 
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bination is produced by first Matte-decomposing the Weyl tensor 

Cabcd = Rabcd + T:R(9ac9db ~ 9ad9cb) 
O 

+ -(gadRcb — QacRdb + QbcRda ~ OdbRca), (3.1) 

into its electric and magnetic parts, 

Eab = C a mbnU m U n , B ab = *C a mbnU m U n , (3.2) 

where U a is a timelike vector field and *C a b c d = \^/~ 9 e abcdC m ( ^, and then 
defining 

T a bcd = C am b n C c j + *C am b n * C c d , (3-3) 
as an energy squared tensor, to give the Poynting vector 

P a = T a bcd U b U c U d 

= C a mbn E mn U b + *C a mbn B m "U b . (3.4) 

Energy tensors of the correct dimension can be constructed using the 
Lanczos tensor [19[. Again there are a large number of possibilities, with 
the added complication of the choice of gauges possible for the Lanczos 
tensor. Here it is assumed that both the Lanczos algebraic 

3 Xa = H a \ = 0, (3.5) 

gauge and differential 

Dab = H a b% (3-6) 

gauge ave been applied. Subject to these gauges the Weyl tensor is recovered 
as a linear diferential equation in the Lanczos tensor 

Cabcd = Habc\d ~ Habd\c + Hcda;b ~ Hcdb;a 

— 9acHb%e + 9adHb.dc;e ~ 9bdH a e C]e + g bc H a e d . e . (3.7) 

The Lanczos tensor has the symmetries 

H a bc Hbaci 

Habc = Hcba ~ H cab . (3.8) 
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Two of the more successfu; energy tensors are 

L a b = H a cd H bcd + *H a c ; d * H bcd 

= 2H a cd H bcd - -g ab H cde H cde , (3.9) 
M ab = 2H a % c . (3.10) 
Their Poynting vectors can be defined as 

O a = L ab U b , 

Q a = M ab U b , (3.11) 

where U a is a timelike vector. 

4 The Metric 

The line element is taken to be of the form 

ds 2 = -dt 2 + dx 2 + dy 2 + dz 2 + h(u,x,y)du 2 , (4.1) 

where 

u = lt + z. (4.2) 

For I = 1, u is null and this is Peres' metric; for I 2 < 1 or I 2 > 1, u is 
spacelike or timelike respectively. For I 2 < 1 define 

Vl - I 2 • z' = It + z = u, Vl-l 2 -t' = t + lz. (4.3) 

Then the metric becomes 

ds 2 = -dt 2 + dx 2 + dy 2 + [1 + (1 - l 2 )h(z', x, y)]dz' 2 , (4.4) 

which is manifestly static. For I 2 > 1 define 

VP - 1 ■ t" = It + z = u, Vl 2 - 1 • z" = t + Iz. (4.5) 

Then the metric becomes 



ds 2 = -[1 + (1 - l 2 )h(Vl 2 - 1 ■ t", x, y)]dt" 2 + dx 2 + dy 2 + dz" 2 . (4.6) 
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This metric has the lapse N as the only varying component of the metric 
and may be written in the form 

ds 2 = -N 2 dt 2 + dx 2 + dy 2 + dz 2 . (4.7) 

When the lapse is of the variables seperable form 

N = H(t)a(x,y), (4.8) 

the time dependence can be absorbed into the metric, showing that the 
metric is static. 

5 The Field Equations 

The line element [4.71 gives the Christoffel symbols 



{^} = f, = §, &y = Nrf, (5.i) 

where i,j = x,y,z and N = d — tN. The Riemann tensor, Ricci tensor and 
Ricci scalar have nonvanishing components: 

R'itj = R i t j t = NN i j , (5.2) 

Ru = NN { { , R ij = -^ (5.3) 
2V-* 

R = -2-^. (5.4) 
The first covariant derivatives of the Riemann tensor are 

Rimt = ^ 2 (^)°> Rm, k = N 2 (^j^ (5.5) 

Rjtik;t = NiN jk -N k Nji. 

The first covariant derivatives of the Ricci tensor can be found by contrac- 
tion; note in particular 

R it;t = N j N ij -N i N j j , (5.6) 
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The second covariant derivatives of the Riemann tensor are 
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itjt;tt 
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+N \m) ~ NN + R W;t) , 
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— (mNjk-NkNji 
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+ N k 



N 



N 



(5.7) 



Rijkl;tt 
Ritjt;kl 



N 



(N — iN — kN 3l - NtNtNjk - N^N^ + NiNjN ik ), 



A7 



The tensors 2.S and |2.6| are 

Z a = 2 (p 



2b 



N 



Ng 



ab 



+ 4(p+- NN*\^ 



N. 



, N 



4bN* 
N 



R. 



it:t 



Zij 



2b 



-N 



R 



it;t 



AT3 J 



4[p+-)N 



N N 



2b 



1 

N 



Ng 



N 



— (NiR^t + NjRu.t] 



(5.8) 



and 

X tt 



p + 



N 



7b 



SbNijN? + I 2p - ^ ) NJN. J 



•i. 'j. 



Vij 



N 



Ng 



ab (N£ 
N 



(5.9) 



(5.10) 
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x it = o, 

respectively. 



6 Solutions to the Field Equations 

From vacuum general relativity, comparing the Ricci tensor |4.3| and the 



Riemann tensor 4.2 we see that the only solution is flat space. For a perfect 
fluid the first and second conservation equations 2.14 become 

M° = 0, 

Gu + ^ + v i = 0, (6.1) 
respectively. The stress is given by 

sH. = -i(3 M + i/), 

Sij = (6.2) 

For the perfect fluid in general relativity, the S xy , S xz , S zy equations show 
that N is a function of one spatial coordinate, say x, and then comparing 
S xx and S yy shows that this vanishes; again the only solution is flat space. 

In general the firld equations constructed in Sec |5| are too complex to 
solve, especially because of the occurrence of Ru-t- Taking 

x = r sin8, y = r cos9, (6-3) 



gives 



N 2 



B, lu = 1 (N N 0>r - ^ - AvA>,, -f 2,-.Y; ) + X,X ri - X r X 



r ly zzi 



Ret-t = N r X e . r - - N e N rr + X z X 9z - X e X zz , (6.4) 

Rzt;t = N r X zr - X z X rr + \[N e N ze - X z (X ee - 2rX r )]. 
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Requiring N = N(t, r) still leaves R r t t t nonvanishing. Ru-t can be made to 
vanish by assuming that 



N = N(t, k — ix % ). 
Then if one lets k 2 = kik 1 the field equations become 



(6.5) 



sl 



A k 2 N" 



K 



k 2 N N 



2b 
\P~ 



Ng ab ( '-tt 



N" 
~N 



+2k z [p + 



N 



Sr 3 



4k 2 K] [p+-)N 



1 

N l~W 



N" 2b 

+ TT 



i 2 N N 



Ng ab ( '-tt 



N 



N 2 



N" 2 + 2NN' 



b 



2 AT"2 



+ 



k z N 



N 2 



,(N_ 

N 



ii\ i 



(6.6) 



N 



+ 



N 



hi 



N 
- —S. j 



N 2 



First consider the case 
b 



p + 



3 6k 2 I m^rp 



rn 



RB , 



and A = 0; then the field equations become 
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2bK 2 



R, 



ab 
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2ZX 



.{U + m 2 RB )R ab . 



RB 



For S a b = 0, by the use of [5^ and |T8| this equation is 







1 /A^ x ° 



N \ N 



N 



N 



m RBNij, 



(6.7) 



(6. 



(6.9) 



where N = d+N. This has the static solution 



N = A(t) exp ( ao + a\Z — -m 2 RB z 2 



(6.10) 
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where ao and a\ are constants. This solution can be rotated around the 
spatial axis, and also transferred into the form 4J. 

A nonstatic solution can be found by assuming that the lapse ./V is a 
function of w = k a x a , where k a is a constant. Symbolically this is 



N = N(-k t + kix 1 ) = N{k a x a ) = N(w). 
The Riemann tensor takes the simple form 



R A 



If one defines 



t = u 



'tj 



It' + z' 



-kk^ 

KiK 3 N ■ 



Z = 



f + y 

the metric can be put in the form JO with 



h = l-N z 



a { —ko + y- ) + />', .(• t kuil + k 



If one assumes that p = —6/3 as in 6.7, the field equations are 



Q t Q I 



2 2 [ 2^" k 2 

f\ TTl^Sij — k*i kj \ TTt 9 7-7- ~ TT 



1 /N'^' 
N \1V 
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+ N 



N 



N N 

When S a b = these equations reduce to the single equation 



1.2 
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11 \ 1 



N \ N 
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+ N 
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N 



ml = 0. 

2 N 



(6.11) 

(6.12) 
(6.13) 

(6.14) 



(6.15) 



(6.16) 



Integrating once, with the constant of integration 01, and multiplying by N 
gives 



(k 2 N 2 - k 2 ) (?pj-rn 2 NN> 



aiN = 0. 



(6.17) 



Setting the constant of integration a\ = and integrating again with the 
constant of integration a 2 gives 



N" 
~N 



1 



a 2 + -m l 2 k~ l \n{k 2 N 2 - fcg). 



(6.18) 
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This differential equation remains intractable when mi ^ 0, but when ni2 = 
it has the simple solution 



N = C + exp(k a x a ) + C- ex.p(-k a x a ) 



(6.19) 



where a factor of ^Ja^ has been absorbed into k a , and thus i/a2~k a — * 
because a,2 is a constant of integration k a may now be complex. For |6.18 
the Riemann tensor takes the simple form 



R' 



itj 



ka k 



J ■ 



(6.20) 



For p + 4 7^ 0, inspection of the S t l term 6.7 and the fluid conservation 



equations If^lj suggests that it should be possible to integrate the field equa- 



tions, but so far this has turned out not to be the case. It is known [20| that a 
sperically symmetric metric is always a solution of the general relativity field 
equations if a sufficiently large number of other fields and fluids are present. 



Adding more terms to the stress in 6.7 has not yet given invertible equations 



To find p + | solutionsthe expresion EL20| , i.e. R 



■.itj 



kok-j 



-N, 



is substituted back into the field equations 6.7. This assumption amounts to 
disregarding higher derivatives of the Rieman tensor and investigating only 
the effect of cross terms; therefore it complements previous linear analysis. 



Substituting |6/7j gives the field equations 



5, 



A: 2 



0. 

kj k^ 



Taking S a b = gives the solution 
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(6.21) 



(6.22) 
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7 Poynting Vectors for the Metric 



From Eqs. |3.1| and 4.2 the Weyl tensor is found to be 



a* 



N- ■ M k 
_LL _l fc- 

' 2N 6N 



1 



(7.1) 



N m 

- 1 'to. 

' 3iV 



The electric and magnetic parts become 



Eij 
B a b 



2N 6N 
0. 



f]iji 



The Poynting vector |3.4| is 

pa 

and thus 



6? 



12N 3 
-1 



p pa _ 



(37V. 2 
(3iV, 



AT-. 2 ) 



2 -n- 2 ) 2 



(7.2) 



(7.3) 



(7.4) 



is always timelike. This Poynting vector appears to give no indication of 
energy transfer as P t is the only nonvanishing component. 

The Lanczos tensor can be calulated by direct methods similar to those 
used in Ref.[21|. In the Lanczos algebraic and differential gauges it is found 
to be 

#. N - ■ 1 



Kit 



- » H m = ^ ( N jVik - N iVjk ). 



The energy tensor 3.1C is 
N? 



Lt 



0, U 



1 



18 ' ^ u v 18iV 2 

and the associated Poynting vector is 

N?6" 



(5NiNj - 2 Vlj N 2 ), 
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18iV 2 ' 



(7.5) 



(7.6) 



(7.7) 
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thus 



0„0" 



\isn 2 J ' 



(7.8) 



and again this is always timelike. 
The energy tensor |3.10| is 
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3 \~N 
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(7.9) 
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3N 2 3N mj 3N ' 



and the associated Poynting vector is 



Q a = — 



~W 



giving 



QaQ a 



~9 9 
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(7.10) 



(7.11) 



For |3g with either C+ = or C_ = 0, Q a = 0; however, for C = C + = C_, 

ft o 

AT2 



QaQ a = ~\k 2 (k 2 - || ) Sec/i 4 ^). (7.12) 



For the solution |6.22|, |7.12| reduces to 



-1 



1 



2A- 



+ ^ ■ ^(^«/). (7.13) 



(™WP _m 2 R B ) \ m WP- m RB 3N2 

When rriRB = 0, Q a is always timelike; when Mrb 7^ a variety of behaviour 
is possible, depending on the values of the constants. 
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